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We study the gluon radiation spectrum off a liard in-medium produced quark in the 
multiple soft rescattering formalism of Baier-Dokshitzer-Mueller-Peigne-Schiff, and 
of Zakharov (BDMPS-Z). Its dependence on the quark and gluon energy, on the gluon 
transverse momentum, on the in-medium pathlength and on the rescattering proper- 
ties of the nuclear medium is analyzed quantitatively. The two components of gluon 
radiation, the hard vacuum radiation associated to the quark production vertex, and 
the medium-induced rescattering contribution interfer destructively. For small spa- 
tial extensions of the medium, this destructive interference overcompensates the hard 
vacuum radiation, and the total medium-induced radiative energy loss decreases as 
IS.E (X — L^. Medium-induced gluon production dominates only above a finite criti- 
cal length L > Lcrit which varies between 3 and more than 6 fm depending on the 
rescattering properties of the medium. Deviations from the BDMPS-L^-behaviour 
persist above Lcrit- The medium dependence of the angular gluon distribution is 
dominated by transverse Brownian k_|_ -broadening. This results in a depletion of the 
low transverse momentum part of both the hard and the medium-induced contribu- 
tion. As a consequence, the medium-induced energy loss outside a finite angular cone 
size can be more than a factor two larger than the total medium-induced radiative 
energy loss. We discuss implications of these results for the jet quenching signal in 
relativistic heavy ion collisions at RHIC and LHC. 
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1. Introduction 



It is a longstanding goal to establish observables which characterize in a model- 
independent way the non-perturbative properties of a nuclear medium. The 
medium-dependence of perturbatively calculable "hard probes" is studied in- 
tensively to this end. A rigorous exarnple available in e-A and p-A collisions is 
the generalized factorization theoremlil of Luo, Qiu and Sterman (LQS) which 
relates the medium-dependence of a class of high-fcf observables to a very 
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small set of higher twist parton distribution functions of the cold target nu- 
cleus. For ultrarelativistic A-A collisions, we are not aware of a corresponding 
factorization theorem for observables which could be measured with existing 
or designed experiments at RHIC or LHC This makes it interesting to search 
outside the framework of higher twist parton distributions for other, rather 
process-independent medium characterizations which allow to distinguish ex- 
perimentally between the non-perturbative properties of hot and cold nuclear 
matter. 

A good candidate for a one-parameter characterization of non-perturbative 
medium properties is the rescattering parameter no C obtained in the QCD 
dipole pictureiH§. It describes the average squared transverse momentum 
transfered from the medium to an ultrarelativistic hard parton within one 
unit of pathlength, tiqC = (q^)/L. This parameter increases with increas- 
ing nuclear density Uq and (as a consequence of the equipartition theorem, 
(q^) oc T) with increasing temperature. However, in order to characterize on 
the basis of uq C differences between hot and cold nuclear matter, one has 
to understand quantitatively the sensitivity of different observables to hq C. 
To this end, we study here the dependence of the medium-induced one gluon 
radiation spectrum off a hard quark on the rescattering parameter no C. 

Medium-induced radiative energy loss and the associated phenomenon of "jet 
quenching" was proposed a decade ago by Gyulassy and Wangii as a sensi- 
tive probe of the non-perturbative properties of excited nuclear matter. Baier, 
Dokshitzer, Mueller, Peigne and Schiff (BDMPS) established!! that its domi- 
nant contribution comes from the rescattering of the emitted gluon and results 
in a quadratic L^-increase of AE/E. An equivalentHi path- integral formula- 
tioni'i was derived independently by Zakharov (Z) - we refer to the formalism 
as BDMPS-Z. Recently, the present author! derived the k^-differential exten- 
sion (|2.1|) of this path-integral formalism which determines the transverse mo- 
mentum distribution of the emitted gluon. Equation (E .11), the starting point 
of the present work, differs from other formulationsE3'lij by a regularization 
prescript ionBB which makes it applicable to finite size media. In particular, 
( |2.1j ) shows in an opacity expansion!! the correct interpolation between totally 



coherent and totally incoherent radiation spectrum for spatially finite media. 
As we recall in section ^, the medium-dependence of is determined by 
the rescattering parameter no C only. 

Quantitatively, little is known about the one gluon radiation spectrum in the 
soft multiple scattering approximation. For the total, k^-integrated energy 
loss AE, numerical results by Zakharovt3 indicate weak deviations from the 
BDMPS-L^-law. However, a systematic study of its dependence on the rescat- 
tering parameter no C and the energy E of the incident quark is missing. 
Moreover, all existing results on AE are obtained by integrating the trans- 
verse gluon momentum up to infinity. In particular, it remains unclear to 
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what extent deviations from the BDMPS-L^-law become more significant if 
the finite kinematical boundary of the k^-integration and the proper regular- 
ization of the gluon radiation spectrum are taken into account. Also, there 
are no results which establish to what extent the main contribution to AE 
stems from the region of small gluon energy ^ = x -C 1 and small transverse 
gluon momentum |kx| -C a; to which the multiple soft rescattering formal- 
ism applies. This motivates a systematic study of the x- and k_|_-differential 
gluon radiation spectrum. The only related quantitative study due to Baier, 
Dokshitzer, Mueller and determines from the k^-differential radiation 

spectrum the phenomenologically important fraction R{Q) of the energy AE 
radiated outside a cone of opening angle 9. However, also in this study, the 
|k_|_|-integration extends to infinity and the questions listed above remain to 
be answered. This is done in the present work which is organized as follows: 

In section 0, we recall the main properties of the one-gluon radiation spectrum 
derived in the soft multiple scattering approximation. Section ^ collects ana- 
lytical expressions obtained in the dipole approximation. Numerical results for 
the gluon radiation spectrum are presented in section ^ and their implications 
are discussed in the Conclusions. 



2. The medium-induced radiation spectrum 

We are studying the medium-induced, k^-differential one gluon radiation spec- 
trum off a hard quark. In the multiple soft rescattering formalism, it reads! 

a(lna;)ak_L [2n)^ Jz- Jyi 

xy c?ue *^^"e ^-^w " ' ' — ■—IC{y = 0,yi;u,yi\uj). (2.1) 

Here, JC denotes the two-dimensional path-integral 

}Cirim),yi;viyi),m\uj) J Prexp rf^ (*|r' " ^^(O^W)] • (2.2) 



The medium-dependence of (|2.1|) stems from the factor n{^) cr(r), where n{C,) 
determines the density of scattering centers in the medium. The dipole cross 
section 

^{r) = 2CA J^2\^o{q±)\' (l-e-'"--) (2.3) 
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contains configuration space information on tlie strengtli of a single elastic 
differential scattering cross section oc |ao(Q±)P- 



The spectrum ( p.lD describes: i) Gluon radiation off a hard nascent quark 
(* = nas) produced during the collision inside the medium at Z- = 0. This 
is the medium-dependence of the QCD analogue of the radiation spectrum 
associated to /5-decay. ii) Gluon radiation off an incoming quark {z^ = — cxd, 
* = in) contained in the wavefunction of the projectile. This is the medium- 
dependence of the Gunion-Bertsch radiation spectrumEl. 



In the configuration-space formulation the quark radiates a gluon at lon- 
gitudinal position i/L {yi) in the amplitude (complex conjugate amplitude). 
The e-regularization of this cross section does not commute with the z± — ±00 
limit which removes the cut-off of the longitudinal integration. To do this an- 
alytically, we consider a medium of finite longitudinal extension positioned 
along the longitudinal axis between and lJ3. We split the longitudinal inte- 
grals into six parts 

z+ z+ L 2+ L L L z+ 2+2+ 

^! 1^1 1^1 1^1 I ■ (2.4) 



UL 2_ z- L VL L L UL 



The six corresponding contributions to the cross section are labelled in an 
obvious way: 

J =^N^C,Y.h- (2.5) 

a[mx) ak_L vr^ jr^ 



Explicit expressions for the /j's, in which the limits z± — > ±00, e — are 
taken, are listed in Appendix ^ Eqs. ( [A.2| )-( [A77| ). They are the starting point 
of the following discussion. We recall several technical points: 

1. Region of validity of (PTT| ) 

Eg. (|2.1| ) is obtained by assuming small angle multiple scattering for both the 
projectile quark and the radiated gluoni. Its region of validity is limited to 
X <^ 1 and |k_|_| <^ cj. For the calculation of the total radiative energy loss, 
integration over x and is required. As discussed below, the x-integration is 
not problematic since the integrand has its main support at small x. The k^- 
integration, however, has significant support in the infrared and the ultraviolet 
regime. 

2. Infrared behaviour of (|2.1|) 

For vanishing medium thickness, the gluon radiation ( |2.1| ) off a nascent quark 
is proportional to and the k^-integration is logarithmically divergent in 
the infrared. To subleading order 0{x), an additional term exp {—iqiyi — Vl)] 
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appears in the integrand of (p.l| ), with q oc the quark massiSi. This term 
provides a natural infrared cut-off. In the following analysis, however, we shall 
not modify by terms subleading in 0{x). We shall rather subtract from 
( |2.1j ) the medium-independent hard radiation associated with the quark pro- 
duction in vacuum. The medium-dependent remainder turns out to be infrared 
finite. This allows us to give a numerical analysis of the medium-induced ra- 
diative energy loss without depending on an infrared cutoff parameter. 

3. Ultraviolet behaviour of ( |2.1| ) 

The kj_-integration over is ultraviolet divergent even if the medium- 

independent part is subtracted off. This is no fundamental problem, since 
the transverse momentum is kinematically bound by the total energy of the 
gluon, |k_L| G [0,a;]. It indicates, however, that the total radiative energy loss 
calculated from ( p.l| ) is sensitive to the high transverse momentum modes 
|k±| = 0{uj) which lie outside the region of validity of its derivation. This is 
discussed in section O. 



3. Dipole approximation 



Two approximation schemes are known for the spectrum ( p. 11) : the opacity ex- 
pansion and the dipole approximation. In the first, one expands the spectrum 
in powers of the density of scattering centersu'BG3li3. This also gives access to 
the radiation off a small finite number of scattering centers. In the dipole ap- 



proximation, one exploits that the integrand of (|2.1|) has its main support at 
small transverse distances r = |r|. This allows to write the dipole cross section 
( |2.3| ) to logarithmic accuracy as 

a{r) = Cr\ (3.1) 



This approximation is used in all discussions of BDMPSH'Hiiy and of Za- 
kharovooEJ'El. It allows to explore the multiple soft rescattering regime in 
which Brownian motion dominates^. Due to ( |3.1| ), the path- integral in ( f^.l| ) 
simplifies to that of a harmonic oscillator 

A 

JCosziy, yhU,yi\fi) = —exp iAB{y^ + u^) -2i Ay -u , (3.2) 
A= ^ B = cos{nAy), (3.3) 

with complex oscillator frequency 
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Here, the dipole strength C appears muhiphed by the homogeneous density no 
of scattering centers, no C. This is the only model parameter in the following 
discussion. 



In Appendix |^, Eqs. ( |A.8| )- (|A.13 ), we write the six contributions Ij to the 
gluon radiation spectrum ( p.5|) in the dipole approximation. For numerical 
evaluation, the integrands are expressed in dimensionless variables by rescaling 
the transverse gluon energy and the longitudinal distance with the modulus 
of the oscillator frequency 



V2\n\ 



V2\n\L. 



(3.5) 



For radiation off a nascent quark, the contributions Ji, /2, /a to (|2.5| ) vanish. 
The remaining three terms can be represented by the diagrams 



14= 



15= 2 Re 



.0 L 



L 



l6= 



L 



(3.6) 



They denote the direct production term I4 for gluon emission inside the target 
in both amplitude and complex conjugate amplitude, the destructive inter- 
ference term I5 which is the characteristic interference term between gluon 
emission in and beyond the target, and the hard vacuum radiation term Iq. 
The analytical expressions read (see Appendix ^ for intermediate steps) 



1 



2ujV2\Q\ 



2Re 



dzi / dz2 exp 



2ki sinh(i±^(5i -^2) 
N 



(2 + 2z)ki cosh(i±^(zi -52; 



N = Z2 sinh 
1 



l + i,^ 



2Re / dz 



ATS 

1 — i) cosh 



—I 



'cosh^ (iiii) 



X exp 



1 + i) k I tanh 



l+i,^ _ 



l + i^ 
z 



(3.7) 
(3.8) 



(3.9) 
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h = ^- (3.10) 



The contribution Jg is exactly (i.e. without taking recourse to the dipole ap- 
proximation) the medium-independent contribution associated to the hard 
radiation off a nascent quark jet propagating in the vacuum!. This follows 
from the leading 0{x) approximation in which only the rescattering of the 
radiated gluon contributes. Subtracting the medium- independent term /g, the 
medium-induced modification to the hard radiation spectrum is given by 1^ 
and /s, 

^^rae, _ ^ ^ h) . (3.11) 



c?(lnx) (ikj_ TT^ 

The radiation spectrum ( |2.5| ) is positive, of course, since one cannot radiate 
less than no gluon. This positivity is guaranteed by the sum J4 + Is + Jq > 
which according to (^) forms a complete square. For the difference ( ^.11 ) 



between the radiation spectrum ( |2.5| ) and the corresponding vacuum contri- 
bution, positivity is not warranted: if less gluons are emitted from the nascent 
off-shell quark in the medium than in the vacuum, then the medium-induced 
deviation turns negative. The interplay between negative [jet enhancement) 
and positive [jet quenching) values of ( |3.11 ) will be studied in what follows. 



4. Numerical Results and Analytical Approximations 

The k^-integration of ( |3.7| )-( p^ is infrared finite. Using polar coordinates 

'0 JO 



Ij = / d(pk k_i_dk±Ij, (4.1) 
Jo Jo 



one finds 



Ia = 47rRe / dzi f dzo \ — 
Jo Jo 1 2 



Jo V21-cosh[(l+i)Az] 

_j_ lg-2Af2sinh[{l+i)A5/2]/2 ^ (^^ 2) 



4 iV2, 

F = 2l-2i coth[(l + i)Az/2] [(1 + i){M^ + Z2) 

+ coth[(l + i)A5/2]] , (4.3) 

'■^ — 1 1 + i 
h = 47rRe / dz 







2 smh[{l + i)z] 
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(1 + ^) e-(l+*)A4^ tanh[(l+i)i/2] \ 

4 smh[(l + i)zl2] cosh[(l + i)z/2] j ' 



Here, x ^ [0, 1] denotes the fraction of the maximal transverse momentum up 
to which we integrate, and denotes the corresponding upper integration 
hmit '^'^ dimensionless rescaled variables 

M, = ^^=. (4.5) 



These expressions determine the x-differential radiation spectrum 

^^med /„ T TP ^ ,\ / J1, med 



-{noC,L,E,x) = 



d{\nx) A^±\<x^ d{\nx) dk± 

= ^N,Cf (h + h) , (4.6) 



which is a function of the rescattering parameter no C, the pathlength L of the 
quark in the medium, the energy of the quark E and the kinematical cut-off x- 
The x-dependence of ( [4.6|) allows to extract k^-differential information from 
the kj_-integrated expression (^]6|), see section [4E| below. The x-integrated 
total radiative energy loss AE is given by 

(noC,L,E,x) = I dxx^^. (4.7) 



E 



A. Analytical Expansion for small L 
1. Radiation off projectile quark 

For the gluon radiation off an incoming quark, a Taylor expansion of the 



radiation spectrum ( p.lD around L = does not exista. One finds 



ds ^ ( f d^^± f 27r , 

NcCf \ I T:-^ ( — \ e 2noCL 



1 + 0(L2)) +0(L2), (4.^ 



rf(lnx)rfkx 7r2 " " (27r)2 KuqC L 
1 (qx-k^)-kx 1 

(q±-kx)^ (qx-k^)2ki ^ki 



which is non- analytic in L. The three contributions in the bracket [] of ( [4.8|) 
combine to 
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r d^q± / 27r \ -^-^rc ^J- 
J (2^0" V^I^J "°"(qx-k^)2ki- ^^-^^ 

This is the Gunion-Bertsch gluon radiation spectrumEl for a target momen- 
tum transfer which is acquired by the transverse Brownian motion of the 
incoming projectile, 

(ql) = 2noCL. (4.10) 



The dipofe approximation does not break down for L where the multiple 
soft rescattering formalism is questionable. Taking L ^ for hqC L = fixed, 
(|4.9| ) matches smoothly onto the single Gunion-Bertsch rescattering result!. 



2. Quark production in the medium 

The medium-dependence ( [4. 6] ) of the gluon radiation spectrum off a nascent 
quark can be expanded in a Taylor series for small in-medium pathlength. We 
find 



h = ^L'Mt + ^L^Ml + 0{L') 



(4.11) 
(4.12) 



Here, the destructive interference term J5 overcompensates for thin media the 
direct production term J4. The resulting medium-induced deviation from the 
jet energy loss in the vacuum is negative, 

h + h = ^L^Ml + 0{L'). (4.13) 



The negative sign of this result indicates jet enhancement, i.e., the medium- 
dependence leads for small L to harder jet fragmentation functions. Shortly 
after the discovery of the EMC effect, Nachtmann and Pirne proposed a 
phenomenological picture which predicts harder fragmentation functions. To 
our knowledge, Eq. ( [4. 13] ) is the first indication for such a jet enhancement 
effect in a QCD-motivated calculation. It will be discussed quantitatively in 
sections 14 Bl-|4 B. 



3. QCD dipole phenomenology 

The rescattering parameter C is the only model-dependent input of the 
medium-dependent gluon radiation spectrum ( p.ll| ). It characterizes the av- 
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erage squared transverse momentum transferee! from the medium to the pro- 
jectile per unit pathlength, see Eq. ( [4.10|) . We study this spectrum for four 
different values of no C: 

1. noC = 0.1 fm"^ = (63MeV)Vfm 

2. noC = 0.5fm-3 = (141MeV)Vfm 

3. noC = l.Ofm-3 = (200MeV)Vfm 

4. noC = 2.0fm-3 = (280MeV)Vfm 

We chose these values after finding in exploratory calculations that they scan 
the region between negligible total energy loss {AE/E < 3 % for no C = 
0.1 fm"^) and very significant energy loss {AE/E > 20 % for hqC = 2.0 fm ). 

Baier, Dokshitzer, Mueller and Schifi^ described the rescattering properties 
of the medium by the transport coefficient q which satisfies q = 2 no C. They 
estimate uqC = 0.05 GeV^ = (500 MeV)^/fm for hot matter at a temperature 
T = 250 MeV. For cold nuclear matter, their estimate uqC = 0.0025 GeV^ = 
(lllMeV)^/fm is based on the relation between the rescattering parameter 
no C and the gluon distribution xG{x). In comparison, the two lower values 
(no C = 0.1 fm~^ and uqC = 0.5 fm~^) of our study scan a reasonable range 
of possible values for cold nuclear matter, and the larger values of n-o C may 
be associated to excited nuclear matter. 

In principle, the parameter n-o C can be determined phenomenologically for 
cold nuc ear matter from the medium-dependence of hard processes. One 
approachi is to relate no C to the parameter Alqs which characterizes in 
the LQS-factorization approachS the amount of transverse momentum trans- 
fered to the hard parton. A first determination from the dijet imbalance re- 
sulted in Alqs = 0.05 — 0.1 GeV^ corresponding to an extremely large value 
UqC = 0.3 — O.6GeV^/fm0. More recent studieJillii use the much smaller 
value Alqs ~ 0-0-'- GeV^. This is still larger than the above estimate based on 
the gluon distribution. Significant experimental and theoretical uncertainties 
remain. Especially, the connection between no C and Alqs assumes a direct 
connection between multiple soft and single secondary (hard) rescatteringll. 
As a consistency check, one should extract no C also from observables sensi- 
tive to the multiple soft rescattering, as e.g. DIS nuclear structure functionsS. 
There, however, existing parametrizations show a non-negligible dependence 
of the dipole strength C on Bjorken Xbj0 or and it becomes important to 
compare the extracted parameters in the corresponding kinematical regimes. 
The present work does not aim at a phenomenological evaluation of existing 
e — A and p ~ A data which takes the above considerations into account. It 
rather scans the parameter space in the region 0.1 fm"^ < noC < 2.0 fm"^ 
which interpolates between cold and hot nuclear matter. 



10 



B. L- dependence of radiative energy loss 



As long as not stated otherwise, we present the k_L-integrated expressions 
( |4.6| ), ( [4 .7] ) for the kinematical boundary x = 1- AH numerical results are for 
a strong coupling constant = 1/3 and A^^ = 3. 
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Fig. 1. (a) Dependence of the medium- induced radiative energy loss ( |4.7| ) on 
the in-medium pathlength L for E = 100 GeV, x = 1 ™d different values of 
the rescattering parameter no C. (b) Double logarithmic presentation of (a) which 
indicates deviations from the BDMPS-L^-law. 



Fig. [^(a) shows the L-dependence of the total radiative energy loss ( [4.7| ) off 
a hard quark of initial energy E = 100 GeV. Irrespective of the value of the 
rescattering parameter tiqC, the total radiative energy loss AE/E increases 
stronger than linear. The double logarithmic plot Fig. 0(b) of the same calcu- 
lation displays deviations from the BDMPS-L^-lawa. 

For sufficiently large L {L < 14 fm), the increase of the radiative energy loss 
is approximately oc L^'^. At smaller values of L (L < 6fm), the radiation 
spectrum shows a numerically very small increase and turns negativ. This is 
consistent with the analytical result ( [4.13|) and prevents us from extending 



the double logarithmic plot Fig. ^(b) into the small L region. 

The sensitivity of the total radiative energy loss on the rescattering parameter 
no C depends strongly on the size of the in-medium pathlength. For L > 8 
fm, a measurement of AE/E with an accuracy of ~ 5% would be sufficiently 
sensitive to distinguish between different values of the rescattering parame- 
ter no C. For smaller in-medium pathlengths, however, this sensitivity is lost 
rapidly. To further analyze the region of negative AE/E and the delayed onset 
of significant total radiative energy loss, we turn now to the x-dependence of 
this spectrum. 
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fraction x fraction x 

Fig. 2. The x-differential gluon radiation spectrum ( [4.6[) as a function of x for 
different in-medium pathlengths L, and different values of the rescattering param- 
eter no C. 

C. X- dependence of the gluon radiation spectrum 

Fig. ^ shows for relatively large in-medium pathlengths L the x-differential 
radiation spectrum ( [4.6| ) which enters the calculation of the total medium- 
induced radiative energy loss (|4.7|) . The spectrum peaks generically at very 
small values of x. This is an important consistency check since the derivation 
of the radiation spectrum (|2.1|) is based on the assumption x <^ 1. 

The maximal radiation probabilities shown in Fig. ^ exceed unity for suffi- 
ciently large pathlength L. This may be taken as an indication that multiple 
gluon emission plays an important role for sufficiently extended targets. 

In Fig. ^, the same x-differential radiation spectrum is plotted for relatively 
small in-medium pathlengths. The distinction between "large" and "small" is 
made empirically by the onset of more complicated structures in the radiation 
spectrum. The only purpose of presenting Fig. ^ separately from Fig. ^ is to 
make these structures better visible. 
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Fig. 3. The same as Fig. |2|, but for relatively small pathlengths L of the quark 
in the medium. The destructive interference between hard and medium-induced ra- 
diation ("jet enhancement") overcompensates the medium-induced gluon radiation 
("jet quenching") for sufficiently small L. 



Irrespective of the value of the rescattering parameter hq C, the x-differential 
medium-induced gluon radiation spectrum Fig. ^ turns negative for very 
small pathlengths L. This is a consequence of the destructive interference be- 
tween hard and medium-induced radiation. As one increases L, more medium- 
induced radiation amplitudes become available. Whether these lead to a 
stronger destructive interference with the hard radiation associated to the 
quark production vertex, or whether they lead to an increasing gluon emis- 
sion probability depends on the quantitative interplay between in-medium 
pathlength L and the rescattering property nQ C of the medium. For hqC = 
0.1 fm~'^ in the upper right panel of Fig. we observe e.g. that the destruc- 
tive interference effect increases up to L < 4 fm. Only for larger pathlengths, 
the medium-induced gluon production probability finally overcompensates the 
destructive interference term and the a;-differential and x-integrated medium- 
induced energy loss turns positive. 

According to Fig. ^, the pathlength at which exact compensation between 
medium-induced radiation and destructive interference occurs shows a strong 
no C-dependence. It decreases from L ^ 5 fm at no C = 0.1 fm~^ via (L ^ 3 
fm, uqC = 0.5 fm"^) and (L ^ 2.5 fm, UqC = 1.0 fm"^) up to (L ^ 2.0 
fm, UqC = 2.0 fm~^). This decrease with increasing uqC is consistent with 
the picture that a significant transverse momentum transfer to the gluon is 
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needed to kick it out of the kinematical region in which destructive interference 
can occur. 

Quantitatively, jet enhancement {AE < 0) at small pathlengths is very small. 
It is almost invisible in the x-integrated radiation spectrum of Fig. |l|. However, 
it delays the onset of a positive jet quenching contribution. As seen from 
Fig. ^, this offset depends significantly on the rescattering parameter tiq C of 
the medium. It will be further quantified in the next subsection. 









= 0.5 fm-3 






rioC = 0.5 fm-3 
L = 1 f m 

— E = 50GeV 


d ln(x) 




d ln(x) 




0.02 


f ' \ 


L 


= 3fm 















0.8 




-- E = 100 GeV 


-0.02 


iiA/ 










E = 200 GeV . 


-0.04 


w \ / 




E = 50 GeV 
E = 100 GeV 


0.4 








\ / 

.. \t \^ 




E = 200 GeV 









0.005 C.C1 0.015 0.02 0.025 0.02 0.04 0.06 0.08 0.1 



fraction x fraction x 

Fig. 4. The energy dependence of the x-differential radiation spectrum results 
in a simple x-scaling ( [4.14 ) 



D. Energy dependence of the gluon radiation spectrum 



In the gluon radiation spectrum the energy E of the initial quark comes 



always multiplied by the fraction x of the energy carried away by the gluon. 
This implies a trivial x-scaling of the radiation spectrum as a function of 
energy, 

7 inas) J (nas) j-, 



(nC, L, x) = -jrr^inC, L, E', -x) . (4.14) 



d{\nx) ' ' ' d{ln.x) ' ' ' E- 



In Fig. ^, we show illustrative examples of this scaling behaviour for a small 
pathlength (where destructive interference effects are important) and for a 
large pathlength. 



If the x-differential radiation spectrum vanishes sufficiently fast at large x, the 
x-scaling ( [4.14| ) results (as a function of 1/E) in a linear increase of the area 
under the x-differential radiation spectrum. Then, the total radiative energy 



loss 



AE 



scales function of energy like 



AE 



E 



{E) = - 
^ ' E 



AE 



E 



{E') 



(4.15) 
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Fig. 5. L-dependence of the fraction of medium- induced total radiative energy 
loss (4.7) for different values of the initial parton energy E and different values of 
the rescattering parameter no C of the nuclear medium. 



As can be seen from Fig. the value of AE/E doubles approximately if 
one scales the initial parton energy by a factor |. This is consistent with the 



scaling law (|4.15|). It holds to good accuracy for all values of the rescattering 



parameter uq C and for all physically realistic pathlengths L: 

Fig. ^ also illustrates a property of the radiation spectrum consistent with 
factorization: One expects on general grounds that for increasing incident 
quark energy and fixed spatial extension of the medium, the medium-induced 
changes to the parton fragmentation become negligible in the limit E ^ oo. 
This is the naive factorization limit in which high momentum spectra depend 
on the initial parton distributions only, and do not test properties of the 
medium. Fig. | shows a significant decrease of the medium-induced fraction of 
radiative energy loss AE/E with increasing jet energy E. However, it indicates 
a kinematical sweetspot where initial parton energies are sufficiently high {E > 
50 — 100 GeV) to make jet quenching a perturbatively calculable quantity and 
yet sufficiently small {E < 1 TeV say) so that the dependence on hq C is 
still significant and makes jet quenching a sensitive probe of non-perturbative 
medium properties. 

Another way of looking at Fig. ^is by plotting the corresponding critical length 
Lcrit above which the total radiative energy loss is larger than one percent. 
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AE/E{L,nt) = 0.01 



(4.16) 
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Fig. 6. The critical in-medium pathlength Lcrit above which the 
medium-induced gluon radiation off the quark exceeds one percent of the total 
quark energy 

As seen in Fig. |], Lcrit increases slowly with the incident energy of the par- 
ton. This is a consequence of the factorization limit: for a medium of finite 
extension, there is an initial parton energy above which medium-induced de- 
viations of the parton fragmentation function become negligible. Also, Lent 
increases obviously if the medium transfers less transverse momentum to the 
rescattering partion. Lent ranges typically between 3 fm and more than 6 fm. 
This is a consequence of the interplay between jet enhancement at small L 
and jet quenching at large L which shifts the onset of a significant radiative 
energy loss to larger values of the in-medium pathlength L. Taking the finite 
geometry of ultrarelativistic heavy ion collisions into accountil, this effect may 
lead to a significant reduction of the jet quenching signal at RHIC and LHC 
since Lcrit is of the same order as the in-medium pathlength of a substantial 
fraction of the produced jets. 

E. Angular contribution to the radiative energy loss radiation spectrum 

The x-dependence of the radiation spectrum ( ^I6D allows to calculate the gluon 
radiation within an angular segment Ak_|_ G [|k_Li|, |k_|_2|], specified by xi = 
(|k_Li|/u;) and X2 = (|k_L2|/u;), 
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Fig. 7. The contributions to the k_|_-integrated radiation spectrum, coming 
from the phase space region |k_L| G [0,xa;] for x = 0.05, 0.15 and 0.25. Arrows 
indicate how the spectrum changes with increasing in-medium pathlength L. 



In Fig. 0, this spectrum is shown for Xi = ^ind small values of X2 = X "C 1. 
For very small transverse momenta |k_|_| < 0.05 cu, we find that the radiation 
spectrum is significantly depleted over the whole range of available gluon en- 
ergies X G [0,1]. This depletion continues to increase with increasing path 
length L even for rather large L, L = 14 fm. Increasing the upper bound on 
the transverse momentum |k_|_| < x^ to x = 0.15 and x = 0.25, one observes 
that the depletion turns into an enhancement for sufficiently large pathlength. 

Even for large x {x > 0.5), a significant depletion of the spectrum occurs 
at small transverse momenta {x = 0.05) but disappears if higher transverse 
momenta are included [x = 0.15). This cannot be the result of destructive 
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interference since it does not show up in the k^-integrated radiation spectrum. 
It comes from the shifting of the gluon transverse phase space distribtion 
which can be understood in terms of a classical muhiple soft rescattering. In a 
classical picture, one expects a characteristic broadening of the hard radiation 

(S), 



as a function of the q^-momentum transfer which grows by Brownian mo- 
tion (q^) oc L. This leads to a depletion of the low momentum region of 
the radiation spectrum without affecting the k^-integrated spectrum. In con- 
trast, the destructive interference between hard and medium-induced radiation 
discussed in section [4 Q and also present in Fig. |^ affects the k^-integrated 
spectrum. 

Comparing the absolute values of the radiation spectrum for x = 0.05, 0.1 and 
0.15 shown in Fig. |^ with the absolute values of the x = 1 radiation spectrum 
in Figs. and |^, one finds that the main contribution to the radiative energy 
loss (|4.7| ) comes from gluons emitted at relatively large transverse momenta 
\kj_\ > 0.25 61;. However, the derivation of the gluon radiation spectrum ( pTip 
assumes |k_L| ^ Jli . We conclude that the calculation of the BDMPS-Z total 
radiative energy loss receives its main contribution from a kinematical region 
in which the validity of the BDMPS-Z approximation is not guaranteed. This 
is not a feature of the dipole approximation studied here: numerical results of 
Gyulassy, Levai and VitevS indicate that also in the opacity expansion the 
dominant contribution to the radiative energy loss is obtained from gluons 
with significant transverse momentum. 



F. Radiative energy loss outside a finite jet 



cone 



From the ^-dependence of the x-integrated radiative energy loss ( [4.7| ), we can 
calculate the quantity 

-^(6 = arcsin(x)) = -^(x = 1) - -^(x) • (4.19) 



This describes the fraction of the medium-induced relative radiative energy 
loss ^ radiated outside a cone of opening angle 9. Our notation in ([4.19|) 
suppresses the dependence of ^(6) on the rescattering parameter UqC, the 
in-medium pathlength L, the quark energy E and the fraction x of the energy 
taken away by the gluon. Numerical results for ( [4.19|) are shown in Fig. ^ for 



relatively large in-medium pathlengths and in Fig. |^ for the small pathlengths 
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for which the destructive interference between hard and medium-induced gluon 
production is significant. 
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Fig. 8. The fraction ( [4.19| ) of the total radiative energy loss AE/E emitted 
outside a jet cone of fixed angle G. 



As seen from Fig. ^ ^(0) does not decrease monotonously with increasing B 
but has a maximum at finite jet opening angle. The reason is that the radiative 
energy loss outside a cone angle G receives additional contributions from the 
Brownian k^-broadening (|4.18|) of the hard vacuum radiation term Jg. Such 



contributions do not affect the total k_L-integrated yield — (© = 0), since they 
result only in a shifting of the transverse momentum phase space distribution 
of the emitted gluon. However, this shift in transverse phase space shows up 
as soon as a finite cone size is chosen. We conclude from Fig. ^ that the total 
k_|_-integrated radiative energy loss^(e = 0) is not the upper bound for the 
radiative energy loss outside a finite jet cone angle -^(9). Depending on the 
rescattering parameter no C and the in-medium pathlength, the latter can be 
larger by more than a factor 2. 



Fig. H displays the angular dependence ([4.19|) for a small in-medium path- 



length where destructive interference becomes important. For the total k^- 
integrated radiative energy loss^(e = 0),thefi gure shows the small negative 
values discussed above and consistent with the analytical result (|4.13| ). Out- 



side a small opening angle B ~ 5°, however, the medium- induced radiative 
energy loss turns positive. This is again the effect of Brownian k^-broadening 
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which populates the high transverse momentum region. 
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Fig. 9. Same as Fig. |8|, but for small in-medium pathlength L = 1 fm. 



5. Conclusion 

In the above analysis of the medium-induced one-gluon radiation spectrum 
(|3.11| ), we found that the hard radiation component associated to the in vac- 
uum production of a hard quark is affected by i) medium-induced radiation 
and ii) medium-induced rescattering. Both effects are non-negligible in a quan- 
titative analysis. 

The first phenomenon is the destructive interference between hard (vacuum) 
and medium-induced radiation. For small in-medium pathlength, analytical 
and numerical results indicate that the hard parton radiates less energy in the 
medium than in the vaccum. This jet enhancement effect is numerically small. 
At larger pathlength, it is overcompensated by a positive medium-induced 
radiative energy loss. However, a quantitatively important remainder of jet 
enhancement is the delayed onset of significant jet quenching. In comparison 
to an unmodified L^-law, the hard parton has to propagate through a medium 
of larger spatial extension L > Lent in order to loose a detectable fraction of its 
total energy. In relativistic heavy ion collisions, the finite geometry limits the 
maximal extension of the in-medium pathlength. The values of Lcrit between 
3 and 6 fermi found in our analysis, may thus imply a significant reduction of 
the observed jet quenching signal. 

The second phenomenon is the Brownian k^-broadening of the hard radia- 
tion due to multiple scattering. This shifts the k^-distribution of the hard 
(vacuum) component to higher transverse momenta. Subtracting the vacuum 
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contribution, one finds that the fraction of the total radiative energy loss R{<d) 
in ( |4.19| ) emitted outside a jet angle B can be more than a factor 2 larger than 



the total medium- induced radiative energy loss obtained for = 0°. In con- 
trast to naive expectation, the total k_L-integrated medium-induced radiative 
energy loss is not an upper bound on the medium-induced radiation outside 
a finite cone. In relativistic heavy ion collisions, this rescattering effect may 
imply a significant enhancement of the jet quenching signal. 

Quantitative statements about the gluon radiation spectrum depend on the 
energy E of the projectile parton, the in-medium pathlength L and the rescat- 
tering parameter no C which is the only model-parameter in our discussion. 
For the spectrum (|3.11| ), we have studied these dependencies in detail. The 



suggested analysis strategy for experimental data is to determine Uq C for 
normal cold nuclear matter phenomenologically and then to compare to the 
values extracted in relativistic heavy ion collisions. The rescattering parame- 
ter can be related in model calculations to the non-perturbative rescattering 
properties of the hot or cold nuclear medium. In this sense, no C is a model 
parameter which can be extracted from experimental data and which allows to 
distinguish between different physical pictures of the excited nuclear medium. 

The present work is restricted to the analysis of the one-gluon radiation spec- 
trum ( p.llj ) derived in the BDMPS-Z formalism. It was one purpose of our 
quantitative study to clarify the kinematical region of validity of this approach. 
In particular, we found that the small-x gluon emission probabilities are very 
high. This indicates the importance of multiple gluon radiation which is not 
contained in the BDMPS-Z formalism. Moreover, the BDMPS-Z formalism is 
based on the assumption of small transverse gluon momentum | k^ | <^ u while 
we find the main contribution to radiative energy loss for |k_|_| = 0(u;). Both 
features question the validity of the BDMPS-Z formalism and await further 
work. We emphasize, however, that the physical origin of the two qualita- 
tively novel effects discussed here, appears to be very generic. Even outside 
the BDMPS-Z formalism, we expect that the destructive interference between 
hard and medium-induced radiation, and the strong implications of Brownian 
kx- broadening for radiation outside a finite jet cone, leave similar, quantita- 
tively important traces in the medium-induced gluon radiation spectrum. 



A. Simplifying the radiation spectrum ( |3.11j ) 



This appendix contains details of how to turn the radiation cross section 
into a numerically accessible form. We consider a nuclear medium of homoge- 
neous density affecting the propagation of the hard quark in the longitudinal 
interval [0,L]. Replacing the in-medium propagator }C{y , yf, u, yi\Lj) outside 
[0, L] by the free propagator /Co, 
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the e-regularization in (|2.1|) can be removed analytically^. Replacing 

the constraint y = in (|2.1|) by a representation of the 5-function, 
/ c?y|^ exp (-ip_L ■ y), one finds 

I, = Re [ [ rfue-m.(kx+px)e-inoa{u) 1 .^2) 
J (27r)^ J p]_ 

xe-|(^-^^^)«o-{u) . ^ /c(r, 0; u, yi\u;) , (A.3) 
p]^ au 

/3 = 2Re I (^^p pVk2^ / e-*P^ '-^-*'^^-'^/C(ri, 0; r2, L\cu) , (A.4) 

^g^ik,.u |_ . i_ ;c(y = 0, 1/,; u, yi\uj) , (A.5) 
(7y au 

/5 = ^y^ ^^Z/L y rfue-^*^-"^-— /C(y = 0,y^;u,LM, (A.6) 

h = X. (A.7) 



In the dipole approximation, these expressions take the form 
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X exp 
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(A.ll) 
(A.12) 
(A.13) 



The variables A, B and D introduced here have different arguments for the 
different terms: 



A2 = A^ 
A,- 
Ai 



2 sin(fi?/L) 



ujVL 



2sin(fiL) ' 



2sin(fi(?/L - yi)) 



D2 = -noC{L - yi) , 



B2 = 


B5 = cos{QyL) , 


(A.14) 


B3 = 


cos(nL) , 


(A.15) 




= cos(fi(?/L , 


(A.16) 




^noC{L - yL) . 
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